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Abstract 

The (ed — > enp) process in the frame of the Bethe-Salpeter approach with a separable 
• kernel of the Nucleon-Nucleon (NN) interaction was considered. This conception keeps the 

covariance of description of the process. Special attention was devoted to a contribution of 
£\j I the L>-states of the deuteron in the cross section of the electrodisintegration. It was shown 

q ' that the spectator particle (neutron) plays an important role. The factorization of a cross 

section of this reaction in the impulse approximation was checked by analytical and numerical 
calculations. 
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1. Introduction 

Study of static and dynamic electromagnetic properties of light nuclei and, espe- 
cially, the deuteron enables more deeply to understand a nature of strong interactions 
and, in particular, nucleon-nucleon interactions. At high energies considerations of a 
nucleus as a nucleon system are not well justified. For this reason the problems to 
study non-nucleonic degrees of freedom (mesons, A-isobars, quark admixtures etc.) in 
an intermediate energies region are widely discussed. However, in spite of it significant 
progress was achieved on this way relativistic effects (which a priori are very important 
at large transfer momenta) are needed to include in the consideration. 
q ■ Other actively discussed problem is the extraction from experiments with light 

nuclei of the information about a structure of bounded nucleons. It requires to take 
into account relativistic kinematics of the reaction and dynamics of interaction. So 
the construction of a covariant approach and detailed analysis of relativistic effects in 
electromagnetic reactions with light nuclei are very important and interesting. 

The electrodisintegration of the deuteron at the threshold has been of interest of an 
investigation for a long time [l]-[6]. The reason is that the electrodisintegration is an 
essential instrument for study a structure of a two-nucleon system. First of all it is an 
electromagnetic structure. The deuteron has been used as a neutron target to get the 
information about neutron electromagnetic form factors. During last 20 years it has 
also been used to receive constraints on available realistic NN potentials. Analyzing 
of the electrodisintegration process we can clarify the role of non-nucleonic degrees of 
freedom. The non-nucleonic effects are often important in few-body systems. And the 
deuteron is one of convenient candidates because complete calculations can in principle 
be performed. 

First experiments were carried out at low transfer momenta but more modern ex- 
periments [7]- [12] performed at high momenta have opened many new questions in a 
region where relativistic effects are important. The experimental results on the differ- 
ential cross section derived from (ed — > enp) reaction are available up to a momentum 
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transfer of about 1 GeV. This situation is very good for investigation of the deuteron 
structure at short distances with the allowance for some exotic effects which have not 
been earlier important. First of all these are the quark degrees of freedom (see [13], 
[14], for instance), but formerly it is necessary to take into account relativistic effects. 

Bethe-Salpeter (BS) approach [15] can give a possibility to consider relativistic 
effects by consistent way [16]. In the paper the deuteron electrodisintegration within the 
covariant BS approach with the separable Graz II interaction kernel is presented. The 
exclusive differential cross section is calculated in the plane wave relativistic impulse 
approximation. 

The paper is organized as follows. In section 2 the relativistic kinematics of the 
reaction and formulae for the cross section are considered. The BS amplitude is pre- 
sented in section 3. The hadron current in the BS formalism is defined in section 4. 
Factorization of the cross section is discussed in section 5. Then the results of our 
numerical calculations are presented in section 6. Finally the discussion of the results 
is performed and further plans are outlined. 

2. Cross section and kinematics 

Let us consider the relativistic kinematics of the exclusive electrodisintegration of 
the deuteron. The initial electron / = (E, I) collides with the deuteron in rest frame 
K = (M d ,0) (M d is a mass of the deuteron). And there are three particles in the 
final state, i.e. electron /' = (E', I') and pair of proton and neutron. In one photon 
approximation (we also neglect the electron mass) a squared momentum of the virtual 
photon q = (u,q) can be expressed via electron scattering angle 9 

q 2 = -Q 2 = (/ - l'f = u0 2 -q 2 = _ 4 |i||i'| sin 2 9 (1) 

np-pair is described by the invariant mass s = P 2 = (p p + p n ) 2 which can be expressed 
through components of photon 4-impulse: 

s = M 2 + 2M d uo + q 2 . (2) 

Lorentz invariant matrix element of the reaction (see Fig. 1) can be written as a product 
of lepton and hadron currents 

M fi = -te 2 (2ir) 4 5 {4 \K - P + q) 

x <l',s' e \f\l,s e > j 2 <np: {P,Sm s )\J^\d: (K,M) >, (3) 

where < /', s' e |j M |/, s e >= u(l' : s' e )^u(l : s e ) is an electromagnetic current (EM). The 
initial (final) electron is described by Dirac spinor u(l,s e ) (u{V , s' e )) . The hadron cur- 
rent < np : (P, Sms)\J fl \d : (K,M) > is a transition matrix element from the initial 
deuteron \d : (K,M) > with total momentum K, projection M to the final np-pair 
\np : (P, Sm s ) > with total momentum P and spin S, projection rris- Thus the unpo- 
larized cross section of the electrodisintegration of the deuteron can be easily written 
as 
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Figure 1: One photon approximation. 



where the factor R connects the final proton angle in the center of mass system (C.M.S.) 
(where the np-pair is rest) with the same in the laboratory system (L.S.): 

V 2 

R= - , — — . (5) 

Vl + r)\p\ - ep^/rj cos dp 

Here p is a momentum of the final proton in L.S., ep = a/p 2 + m 2 , 6p is an angle 
between the final proton and Z-axis, m is a nucleon mass, and r\ = q 2 /s. 
Tensor of unpolarized leptons in (4) is expressed as 

l ^ = \j2 < r ' £\3l*% S e >< I, S e \ju\l', S'e >= 2(^Z„ + W + 9^ ( 6 ) 

S e S e 

and hadron tensor can be written as 

w p» = \Y< <d:(K,M)\J^\n P :(P,Sms)><n P :(P,Sm s )W-(K,M)> . (7) 

MSrris 

In order to average on initial and sum on final states let us introduce a helicity tensor 
which can be directly connected with structure functions (see, for example [1],[4],[6]). 
These quantities allow to calculate polarization and asymmetry observables easily and 
may be necessary in future (we didn't calculate them in this work). Keeping in mind 
the Hermitian properties of the lepton and the hadron tensors the cross section can be 
rewritten as 

d b <J Cr Mott j-j-, 



dE'dVt'dVtp 8M d (2vr) 3 
x [l° 00 W 00 + l° ++ (W ++ + W__) + 1% _2ReW+_ - l° 0+ 2Re(W 0+ - W .)] , (8) 



where a Mo tt = {ot cos %/2E sin 2 |) 2 is Mott cross section for point-like particles and 



,o _Q 2 ,o Q Q 2 ,^ ,o K^o.Q 2 ,o _Q 2 

' 00 "g 2 ' lo+ -\q\V2Vq 2 +tg 2' '++ ~ 2 %g 2 + 4g 2 ' l+ ~ ~ 2q 2 ' 
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So the calculation of the cross section (8) comes to the calculation of the hadron tensor 
W^u which describes the NN interaction and is a subject of our investigation. 

3. Bethe-Salpeter Amplitude of the Deuteron 

In the Bethe-Salpeter approach (BSA) the deuteron is a bound system which can 
be described by the amplitude $m(&; K) of the equation 

<W(£;; K) = iS§ + kj Sf p (y - *) / ^*Wa(*> ^ ^)<W^; K), (10) 

here S^(K/2 — (— l)^fc) is a propagator of the £-th nucleon, V(k, k'; K) is a ker- 
nel of a NN interaction (Greek letters means spinor indexes). The amplitude of 
the deuteron in rest frame can be expanded through two-nucleon relativistic states 
\aM) = |tt, 23+1 L P jM): 

$ M (k;K) = J2Mko,\k\)y aM (k), (11) 

a 

where S 1 denotes a total spin of the system, L is an orbital angular momentum, J is 
a total angular momentum with a projection M. Quantum number p counts positive- 
and negative-energy states, it marks the parity of the state; (p a (k , \k\) is the radial 
part of the BS amplitude. The spin-angular part y a M(k)) is 

y a M(k)u c = 

* L E ^ZfrC^SrnsC^Y^ (fc)< (p)< T (-fc). (12) 

m L msmim 2 pip2 

Using partial wave decomposition (11) we can write the decomposed Bethe-Salpeter 
equation for the radial part of the BSA 

Mko,\k\) = S a (k ,\k\;s) I dk' j k' 2 d\k'\J2V ab (k ,\k\,k' ,\k'\;s)MKAk'\). (13) 

J J b 

Then taking into account separable anzats for the rank kernel of interaction 

N 

V ab (k , \k\X, \k'\; s) = X v 9?(k , \k\) g b j(k' , \k'\), X tJ = \ jh (14) 
we can find a solution of the BS equation (13) in the following form 

N 

<j> a {ko,\k\) = ^2 S a (k ,\k\; s) Xij gi(k ,\k\) Cj(s). (15) 

Here \j are fitting parameters, gf(ko, \k\) are trial functions. The coefficients Cj(s) 
satisfy the following system of homogeneous equations 

N 

c i( s ) ~ E h ik( s ) ^kj Cj(s) = 0, (16) 
k,j=l 
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where hik(s) are defined by the integral 



h ^ s ) = ^H J dk ° I k 2 d\k\S a (k ,\k\;s)gt(k ,\k\)g a k (k ,\k\). (17) 

a 

Using the covariant separable Graz II rank III kernel of interaction [16] we can find \j 
(see Table 1) from an analysis of the experimental data for the deuteron characteristics 
(phase shifts, binding energy, length of scattering etc.). 



Table 1: Parameters of the covariant separable Graz II kernel 
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It is necessary to remark that here we took into account the positive-energy states 
( S S+ s £>+) only 



9 ?Hko,\k\) = ;- ( i 8) 
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g 2 ^ (k , \k\) = - /> 9 VU , 2 ^ (19) 



gl Dt (k , \k\) = gl Dt (k , \k\) = gl St (ko, \k\) = 0. (21) 

In our calculation it is more convenient to use the BS vertex Ym which is related 
with the BS amplitude by simple expression. For full functions 

$ M (k; K) = fifW (| + kj S< 2 > - kj T M (k; K), (22) 

and for their radial parts 

<p a {k , \k\) = ^ Sab{k , \k\;s)g b {k , \k\), (23) 

b 

here we omitted spinor indexes for simplicity. S a b is a nucleon propagator which is 
diagonal for positive-energy partial parts 

S ++ (k , |fc|; s) = 1/ [(v^/2 + k - e k ) (^i/2 - k - e fe )] . (24) 
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Figure 2: Relativistic impulse approximation. 



4. Hadron Electromagnetic Current 

Let us write the matrix element of the hadron electromagnetic current with the BS 
amplitude using Mandelstam technique [17] 

< np : (P, Sm s )\J ll \d : (K, M) >= 

J0f (Si (p; v\ P) Mp, fc; p, k) $ M (k; K). (25) 

Here p* is the relative 4-momentum of the np-pair (nucleons are on-mass-shell), Pp* = 
0, and p* 2 = s/4 - m 2 . 

We consider the process of the electrodisintegration of the deuteron in RIA (see 
Fig. 2). In our further calculation only one-body currents are taken into account 

AW (p, k; P, K) = *(2vr) 4 j*<*>(p - k - |)r« + p, | + *) 5^ - p) 



p 



p 



■(26) 



(here the total and relative momenta are introduced: P = p\ + p%, K = ki + k2, k 
l(ki - k 2 ),p = i(pi - pa))- 

In this case the matrix element of the hadron current has the following form 

d 4 p 



<np:(P,Sm s )\J^\d:{K,M) >= i J] / 

£=1,2 ^ 



if (g) 



P 



(2?r) 4 X5m s (p;P) 

r M (p + (-i)^;^). 

Note that jNN- vertex was taken on-mass-shell 

r<fV,p) — rW(g) = 7 ,pf V) - ^ ( 7 ,g - qi») pf(g 2 ). 



(27) 



:28) 



Here (F% f') - Dirac (Pauli) form factor of the nucleon which obeys the next nor- 
malization conditions 



if'(0) 



F 2 w (0) = *i 



■p; 



P 1 (2) (0) = 0, F 2 w (0) = xi 



(2), 



(29) 
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Figure 3: Plane wave approximation. 



(><n) is the anomalous proton (neutron) magnetic moment. 

It should be mentioned that we neglect a final state interaction (FSI) (it is a subject 
of future calculations). 

= (2-) 4 * W (P - Pi E C l7 lh m 2 ^ (f + ^ (f - P) ■ (3°) 

mirri2 

Using (30) for the final np-pair wave function and integrating (27) over p we obtain 
our basic RIA expression for the hadron current 

< up : (P,Sm s )\JM : (K, M) >= i £ x£ ff (p*,P)rj?( ff ) 

£=1,2 

xs w (f -p* - Tm ( p - + (-i)'f ; *) • (31) 

It has very simple form. And to get it we should just perform the analytical calculation 
of the trace. For this purpose we use the REDUCE system. 

5. Factorization of the cross section 

Let us consider the electrodisintegration of the deuteron supposing that an initial 
lepton collides only with the proton in the deuteron and the neutron is a spectator. 
Then the cross section is factorized on two parts, one is connected with the contribution 
of the neutron as a spectator and another with a proton contribution, the latter does 
not have interference terms between the S- and D-states. 

5.1 Nonrelativistic case 

The amplitude of the process can formally be presented as a production 

M=xt 1 xl 2 0^ M , (32) 

where spinors xti x > Xm 2 describe the outgoing np-pair, O corresponds to the interaction 
vertex, $m is a wave function of the deuteron. Let us note that the vertex O stands 
in general for any one-particle interaction, but in this paper describes 7iViV-vertex. 
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Inserting into this expression a complete set of pair states we can get 



(33) 



after evident transformations the hadron tensor can be written as 
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Introducing the partial-wave decomposition for the deuteron: 

lms\S2S 

and transforming the second term with the help of orthogonalization properties of the 
spinor \ an d some relations for Clebsh-Gordan coefficients we can finally obtain the 
factorized expression 



W, 



^K\Y\ui\\ 



(34) 



where = xtiflnXm^ is a one-body interaction part and / counts partial states of the 
deuteron. 

Thus it is seen that the cross section is proportional to the sum of squared radial 
parts of the deuteron wave function. 

5.2 Relativistic case 

In the relativistic case the matrix element of the deuteron electrodisintegration can 
be written schematically in the following form 



M = 



pair 



o®s®r M , 



(35) 



where ^ pa i r is the wave function of the np-pair, O is a vertex of interaction, S is a 
propagator of the nucleon, Ym is the vertex function of the deuteron. 

Introducing the partial-wave decomposition of the deuteron vertex function in the 
L.S., considering the only proton interacting with a virtual photon and supposing the 
PWA for the final np-pair we can present the comprehensive expression in the following 
form 

s\S2 m'[ 

x Y C l™ m Cr^ m u {1 \rni, ~kM 2 \m2, -k 2 )Y lm (k) 9l (k 0: |fc|), 



m\m2lmsm a 



where S\ + {ki) = 1 / (fcio — e^, ) and T^\q) vertex is described by Eq. (28). As it was 
assumed above only 3 S±, 3 -D^-states were taken into account. Using orthogonalization 
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properties of the bispinors and some relations for Clebsh-Gordan coefficients we can 
write 

s\S2m\lmsms 

with 

4 1 )(si,p 1 ;mi,fei) = ««(si,p 1 )rW(g)««(mi,fci), (36) 
is a one-body photon-proton interaction part. Now we can derive the hadron tensor 

w v> = \ E M,Mu. 

MSM S 

Using once more properties of the Clebsh-Gordan coefficients and Dirac spinors we 
obtain the expression 



W, 



flU 



C d Sp^p, + m)T^(q)(ki + m)f W(q)} (37) 



which involves the simply calculated trace and a function 

C d = -L-J—S^fcO \9i(ko, \k\)\ 2 

containing the structure of the deuteron. Performing the trace calculation we finally 
obtain the expression for the hadron tensor 

= c d {w; u F?( q 2 ) + w%H?W) + (38) 

with 

= 4 [Pi^iu + k lfM p lu + [m 2 - ( Pl ■ h))g^] (39) 
W* u = 2 [ki^q u - q^k lu - pi^q u + q^pi u ] 

= [(-q 2 m 2 - q 2 ( Pl ■ k x ) + 2( Pl ■ q)(k 1 ■ q)) + (m 2 + ( Pl ■ h)) q^q v 
~ {(ki ■ q)(pi^qu + %Piv) + (Pi • q)(ki^q u + q^k lv ) - q 2 (pi fl k ll/ + k ltJl p lv ))] /m 2 . 

Let us note here in the expressions Eqs. (37,39) the four-vector k\ has the on-mass-shell 
form k\ = (e^ , ki) in differ with k\ = (k w , ki) in the Fig. 2. 

Thus we see that the factorization of the electrodisintegration cross section exists 
both in nonrelativistic and relativistic cases. The necessary conditions for this are 
the plane-wave approximation for the final np-pair, the neutron in the deuteron is 
supposed to be a spectator (the one-body type of the interaction in the vertex O) and 
only positive-energy states for the deuteron are taking into account. As for the second 
condition the type of one-body interaction does not play any role but only spin-one-half 
particle is scattered. The third condition means the P waves in the deuteron (namely 
3 P^~ and 1 Pi~) destroy the factorization. 

6. Results and Discussion 

We present here the results of the calculation of the deuteron electrodisintegration 
cross section in the relativistic plane wave impulse approximation with the separable 
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Figure 4: The electrodisintegration cross section versus the neutron momentum for three 
kinematics of the experiments at Saclay. Solid and dashed lines correspond to the calculations 
with and without neutron contribution (upper plot). Bottom plot shows the relative neutron 
contribution in the corresponding experimental regions. The experimental data regions were 
taken from [7](Saclayjji) and [8](Saclaym). 

Graz II rank III kernel of interaction. In our calculations we follow to conditions of 
real experiments and we distinguish eight sets of experimental data. Let us mark these 
sets as Saclay i, Saclay n (see [7], Table 3); Saclay ui ([8], Table 1); Bonnie Banna 
(see [9], Table 3); Bonnjjj, Bonniv, Bonny (see [12], Tables 5,3,4 respectively). 

First of all we illustrated the influence of the spectator neutron on the cross section 
(see Figs. 4, 5, 6). It is seen that it increases with the increasing of the neutron 
momentum and reaches 50% in the Saclaym kinematic range. One can see that 
the cross section of the deuteron electrodisintegration versus \fs changes not so strong, 
nevertheless the contribution of the spectator neutron is not negligible. Let us note that 
this contribution changes sign in the Saclaym kinematic region (see Fig. 4). In order 
to understand the origin of this behavior we present on the Fig. 7 partial contributions 
of the S- and D-states for this kinematical region versus neutron momenta. We found 
that the D-state plays an important role and then it is naturally to ask what happens 
if we change the magnitude of the D-state. On the Figs. 8, 9, 10 we can see that 
for different magnitudes of the D-states the cross section changes distinctly (especially 
for the kinematic region [8]), but the ratio 5 = (a p+n — a p )/a p+n is not changed at 
all (see Fig. 11). It means that the difference is mainly connected with the spectator 
neutron contribution. Thus we can make a conclusion that the experimental data 
within the kinematics from Saclaym [8] can supply the good test for various models 
of NN interactions in the deuteron. 
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Figure 5: The same as in previous figure but versus pair invariant mass yfs for the kinematical 
conditions were taken from [9] (Bonnjjj). 




Figure 6: The same as in previous figure. The kinematical conditions were taken from 
[12}(Bonn in j V y). 
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Figure 7: The contributions of the spectator neutron versus the outgoing neutron momentum 
to the electrodisintegration cross section of the deuteron partial S-, D-states are shown for 
the experimental sets from [8](Saclayiu). 




Figure 8: The contributions of the deuteron partial D-state to the electrodisintegration cross 
section versus neutron momenta are shown for three sets of the experiments [7](Saclayiji), 
[8] (Saclayin). 
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V s, GeV 

Figure 9: The contributions of the deuteron partial D-state to the electrodisintegration 
cross section versus pair invariant mass y/s are shown for the conditions of the experiments 
[9] (Bonn/,//). 




V s, GeV 

Figure 10: The same as in the previous figure but for [12] conditions (Bonnjuy). We 
omitted curves for the Bonniv case because they are very close to Bonny- 



13 




n (MeV/c) 



Figure 11: The contribution of the spectator neutron versus neutron momenta to the electro- 
disintegration cross section for different deuteron D-states for conditions of the experiment 
Saclayni [8]. In the first picture solid (dashed) line stands for p + n- (p-) contribution with 
different D-states in the deuteron: po = 4% for lower line and pn = 6% for upper line. 

7. Summary 

In the presented paper we have considered the electrodisintegration of the deuteron 
in the Bethe-Salpeter approach. It is realized for the two-nucleon system by using the 
multipole expansion with the spinor structure of two nucleons. The separable ansatz 
for the interaction kernel has provided a manageable system of linear homogeneous 
equations for deriving the BS amplitude. 

We have switched then to the case with the use of the covariant revision of the Graz 
II separable potential with the summation of several separable functions. 

The reaction of the deuteron electrodisintegration served as a testing ground for 
the method under investigation and helped to outline both strong and weak points of 
the approach. The analysis has proved the technique to be very promising, even if 
we find an evident discrepancies with experimental data at this stage of development. 
Several items can be suggested for the program of further theoretical study. First 
of all it is necessary to take into account the final state interaction for the np-pair. 
Then we need to take into account the negative-energy states for the BS amplitude 
and calculate the contribution of the P waves in the electrodisintegration. After that 
we can calculate different asymmetries of the (ed — > enp) process which can give new 
qualitative information about the structure of the deuteron. 
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